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An overview of the concept of phase-locking at the non linear, geometric and quantum 
level is attempted, in relation to finite resolution measurements in a communication 
receiver and its 1// noise. Sine functions, automorphic functions and cyclotomic arith- 
metic are respectively used as the relevant trigonometric tools. The common point of 
the three topics is found to be the Mangoldt function of prime number theory as the 
generator of low frequency noise in the coupling coefficient, the scattering coefficient and 
in quantum critical statistical states. Huyghens coupled pendulums, the Adler equa- 
tion, the Arnold map, continued fraction expansions, discrete Mobius transformations, 
Ford circles, coherent and squeezed phase states, Ramanujan sums, the Riemann zeta 
function and Bost and Connes KMS states are some but a few concepts which are used 
synchronously in the paper. 
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1. Introduction 

The interleaving of frequencies and phases of electronic oscillators interacting in 
non linear circuits follows arithmetical rules. Continued fraction expansions, prime 
number decompositions and related number theoretical concepts were successfully 
used to account for the experimental effects in mixers and phase- locked loops |^3I2I ■ 
We also made use of these tools within the field of quantum optics emphasizing 
the hidden connection between phase-locking and cyclotomy [2]. As a matter of 
fact the understanding of phase effects in devices is so much entangled that it will 
reveal useful to check several clues relying on differential and discrete non linear 
equations, Fourier analysis, hyperbolic geometry and quantum mechanics. But 
over all the fundamental issue is the modelling of finite resolution measurements 
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and how it puts a constraint on the performance, how it generates the fluctuations 
in oscillating circuits. 

In Sect. |21 we report on the early history of phase- locking about classical ob- 
servations of coupled mechanical pendulums and electronic oscillators. Non linear 
continuous and discrete generic models are introduced emphasizing the case of ho- 
modyne detection in a phase sensitive communication receiver. It is shown that 
low frequency noise happens due to low pass filtering and the finite resolution of 
frequency counts. A phenomenological model relating the coupling coefficient to 
prime numbers is developed. In Sec the justification of the model is given by 
using the hyperbolic geometry of the half-plane applied to the low pass filtering. In 
Sec 0] the phase itself is taken to be discrete. Discrete quantum optics is related to 
quantum phase-locking and the arithmetic of 1// noise. 

2. Classical Phase-Locking: from Huyghens to the Prime Numbers 

Being obliged to stay in my room for several days and also occupied in making 
observations on my two newly made clocks, I have noticed a remarkable effect which 
no one could have ever thought of. It is that these two clocks hanging next to one 
another separated by one or two feet keep an agreement so exact that the pendulums 
invariably oscillate together without variation. After admiring this for a while, I 
finally figured out that it occurs through a kind of sympathy: mixing up the swings 
of the pendulums, I have found that within a half hour always return to consonance 
and remain so constantly afterwards as long as I let them go. I then separated them, 
hanging one at the end of the room and the other fifteen feet away, and noticed that 
in a day there was five seconds difference between them. Consequently, their earlier 
agreement must in my opinion have been caused by an imperceptible agitation of the 
air produced by the motion of the pendulums. 

The citation is taken from [I]. The authors remind a later letter by Huyghens 
that the coupling mechanism was in fact a small vibration transmitted through the 
wall, and not movement of air: 

Lord Rayleigh (1907) made similar observations about two driven tuning forks 
coupled by vibrations transmitted through the table on which both forks sat. . . Locking 
in triode circuits was explained by Van der Pol (1927) who included in the equation 
for the triode oscillator an external electromotive force as given in 

d 2 v d 

where g is the linear net gain (i.e. the gain in excess of losses, f3' the saturation 
coefficient, and uj is the resonance frequency in the absence of dissipation or gain. 
He showed that when an external electromotive force is included, of frequency loq, 
and tuned close to the oscillator frequency uj, the oscillator suddenly jumped to 
the external frequency. It is important to note that the beat note between the two 
frequencies vanishes not because the two frequencies vanish, not because the triode 
stops oscillating, but because it oscillates at the external frequency. 

We can show the locking effect by utilizing the slowly varying amplitude approach, 
including a slowly varying phase $ and oscillation at the external frequency loq and 



amplitude V 

— — + ifsin<I> = uj — ujq — lolf, (2) 
where we use ojlf for the detuning term and K — loqVq/V for the locking coefficient 

The regime just described is the so-called injection locking regime, also found 
in injection-locked lasers. The equation (j2J is the so-called Adler's equation of 
electronics jH]- 

One way to synthesize © is thanks to the phase-locked loop of a communication 
receiver. The receiver is designed to compare the information carrying external os- 
cillator (RF) to a local oscillator (LO) of about the same high frequency through a 
non linear mixing element. For narrow band demodulation one uses a discriminator 
of which the role is first to differentiate the signal, that is to convert frequency mod- 
ulation (FM) to amplitude modulation (AM) and second to detect its low frequency 
envelope: this is called baseband filtering. For more general FM demodulation one 
uses a low pass filter instead of the discriminator to remove the high frequency sig- 
nals generated after the mixer. In the closed loop operation a voltage controlled LO 
(or VCO) is used to track the frequency of the RF. Phase modulation is frequently 
used for digital signals because low bit error rates can be obtained despite poor 
signal to noise ratio in comparison to frequency modulation 6 . 

Let us consider a type of receiver which consists in a mixer, in the form of a 
balanced Schottky diode bridge and a low pass filter. If /o and / are the frequencies 
of the RF and the LO, and 0{t) and ip(t) their respective phases, the set mixer and 
filter essentially behaves as a phase detector of sensitivity uq (in Volts/rad.), that 
is the instantaneous voltage at the output is the sine of the phase difference at the 
inputs 

u(t) =n o sin(0(i) - ip(t)). (3) 

The non linear dynamics of the set-up in the closed loop configuration is well de- 
scribed by introducing the phase difference <I>(i) = sin(8(t) — if>(t)). Using 9 = ujq 
and ip(t) = uj + Au(t), with ujq = 2irfo, uj = 2irf and A (in rad. Hz/Volt) as the 
sensitivity of the VCO, one recovers Adler's equation @ with the open loop gain 
K = u a A. 

Equation J2J) is integrable but its solution looks complex [J]. If the frequency 
shift ujlf does not exceed the open loop gain K, the average frequency (4>) vanishes 
after a finite time and reaches the stable steady state $(oo) = 2ln + sin" 1 (lulf/ K ), 
I integer. In this phase-tracking range of with 2K the RF and the LO oscillators 
are also frequency-locked. Outside the mode-locking zone there is a sech shape beat 
signal of frequency 

Q LF = {$>{t)) = (Lo 2 LF -K 2 ) 1 ' 2 . (4) 

The sech shape signal and the non linear dependance on parameters lolf and K are 
actually found in experiments [HE]- In addition the frequency ujlf is fluctuating 
(see Fig. 1). It can be characterized by the Allan variance <t 2 (t) which is the mean 
squared value of the relative frequency deviation between adjacent samples in the 
time series, averaged over an integration time r. Close to the phase-locked zone the 
Allan deviation is 

°{r) = (5) 

^LF 
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Fig 1. Fluctuating counts of the beat frequency (in Hz) close to the phase-locked zone. The inputs 
are quartz oscillators at 10 MHz. The power spectrum has a pure 1/f dependance 

where cxo is a residual frequency deviation depending of the quality of input oscilla- 
tors and that of the phase detector. Allan deviation is found independent of r which 
is a signature of a 1/f frequency noise of power spectral density S(f) = c/(2 In 2/) . 
One way to predict the dependence JSJ is to use differentiation of J3J with respect 
to the frequency shift ujlf so that 

8Cj lf = 5lo lf {1 + KyCulp) 1 / 2 . (6) 

Relation © is defined outside the mode-locked zone \u>lf\ > K; close to it, if the 
effective beat note ujlf < K, the square root term is about K/ujlf. If one identifies 
5ujlf/&lf with a pedestal Allan deviation <jq and SQlf/ijJlf with a magnified Allan 
deviation a one explains the experimental result (JjjJ. One can conclude that, either 
the PLL set-up behaves as a microscope of an underlying flicker floor ao, or the 1// 
noise is some dynamical property of the PLL. In the past we looked at a possible low 
dimensional structure of the time series and found a stable embedding dimension 
lower or equal to 4 But at that time the dynamical model of 1// noise still 
remained elusive. 

Adler's model presupposes a fundamental interaction ujlf = |k>o — m thc- 
mixing of the two input oscillators. But the practical operation of the phase detector 
involves harmonic interactions of the form ujlf — \p^o — < uj c — 2-7r/ c , where 

p and q are integers and f c is the cut-off frequency of the low pass filter. This 
can be rewritten by introducing the frequency ratios v = and /i = as 
fi = q\v — ^|. This form suggests that the aim of the receiver is to select such 
couples (p, q) which realize a "good" approximation of the "real" number v. There 
is a mathematical concept which precisely does that: the diophantinc approximator. 
It selects such couples pi and qi, coprime to each other, i.e. with greatest common 
divisor (pi,qi) = 1 from the continued fraction expansion of v 

Pi 

v = [a ;ai,a 2 , ■■■a i ]=a + l/(ai + l/(a 2 + 1/ h l/(a; ■ ■ •))) ^ — ■ (7) 



The diophantinc approximation satisfies 



W- P -\<^. (8) 

Qi a i+ iqf 

The fraction ^ is a so-called convergent and the a'^s are called partial quotients. 
The approximation is truncated at the index i just before the partial quotient aj+i. 
It should be observed that diophantine approximations arc different from decimal 
approximations §*■ for which one gets \v — < -r. It was shown using the 
filtering condition that aj+i identifies with a very simple expression 



a-i+i 



Jo_ 



(9) 



where [ ] denotes the integer part. For example if one chooses /o = 10 MHz and 
f c = 300 kHz, the fundamental basin ^ = 1 will be truncated if a; + i > 33 and 
the basin jjr = § will be truncated if aj+i > 6. The resulting full spectrum is a 
superposition of V-shape basins of which the edges are located at 

V\ = {°0; 7 Oli a 2, ■ ■ ■ 7 Oi, Oi+i}, 

^2 = {a ; ax, a 2 , • • • , Oi_i, 1, a^+i}, (10) 

where the partial expansion before a^+i corresponds to the two possible continued 
fractions of the rational number 2i. The basin of number v = | = {0; 1, 1, 2} ex- 
tends to z/i = {0; 1, 1, 2, 33} = |§ ~ 0.594, i/ 2 = {0; 1, 1, 1, 1, 33} = § ~ 0.618. For 
a reference oscillator with /o = 10 MHz this corresponds to a frequency bandwidth 
(0.618 - 0.594). 10 7 MHZ=240 kHz. 

With these arithmetical rules in mind one can now tackle the difficult task to 
account for phase- locking of the whole set of harmonics. The differential equation 
for the phase shift Q(t;qi,pi) at the harmonic (pi,qi) corresponding to the beat 
frequency 

u L f = \piUo - g»w(f)|) (ii) 

can be obtained as 

^(t;q l ,p l )+q l H(P) ^K(r uSl ) 

xsin(— $(t;qi,pi) - — (cfcrj - PiS t ) + $o(«"i,Si)) = W L F(p»,?i)< (12) 
qi qi 

The notation K(ri,Si) means the effective gain at the harmonic rj/sj, ( E'o( r i7 Si) is 
the reference angle and H(P), where the operator P = 4i, is the open loop transfer 
function. Solving l|12l) is a difficult task. Let us observe that the RF signal at 
frequency ojq acts as a periodic perturbation of the Adler's model of the PLL. If 
one neglects harmonic interactions, Ijl2(l may be simplified to the standard Arnold 
map model 

$„+i = $„ + 2nn - c sin (13) 

where Q = j-j- is the bare frequency ratio and c = -J^. Such a nonlinear map is 
studied by introducing the winding number v = lim„_, 00 ($ n — $o)/(27r^)- The limit 
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exists everywhere as long as c < 1, the curve v versus w is a devil's staircase with 
steps attached to rational values Q = ^ and width increasing with the coupling 
coefficient c. The phase-locking zones may overlap if c > 1 leading to chaos from 
quasi-periodicity [§]• 

The Arnold map is also a relevant model of a short Josephson junction shunted 
by a strong resistance R and driven by a periodic current of frequency uq and 
amplitude I - Steps are found at the driving voltages V r — RIq — r(hiu /2e), r 
a rational number. Fundamental resonances r = n, n integer, have been used to 
achieve a voltage standard of relative uncertainty 10~ 7 . 

To appreciate the impact of harmonics on the coupling coefficient one may ob- 
serve that each harmonic of denominator qi creates the same noise contribution 
Sujlf — qiSui(t). They are 4>{qi) of them, where <f>(qi) is the Euler totient func- 
tion, that is the number of integers less or equal to qi and prime to it; the average 
coupling coefficient is thus expected to be l/0(q<j). In [2] a more refined model 
is developed based on the properties of prime numbers. It is based on defining a 
coupling coefficient as c* = cA(n; qi,pt) with A(n;qi,pi) a generalized Mangoldt 
function. Mangoldt function is defined as 



The generalized Mangoldt function attached to the resonance 2i adds the restriction 
that n should also be congruent to pimod(q;). The important result of that analysis 
is to exhibit a fluctuating average coefficient as follows 



with e(t) = 0(t~ x l 2 hi 2 (i)) which is known to be a good estimate as long as < y/i 
PQ. The average coupling coefficient shows the expected dependance on q^. In 
addition there is an arithmetical noise e(t) with a low frequency dependance of the 
power spectrum reminding 1/f noise. Although that stage of the theory is not the 
last word of the story, it is quite satisfactory that this approach, based on phase- 
locking of the full set of harmonics, is accounting for the main aspects of 1// noise 
found in experiments. 

3. Hyperbolic Phase-Locking 

The concept of an automorphic function is natural generalization of that of a peri- 
odic function. Furthermore, an automorphic form is a generalisation of the expo- 
nential function e(z) = exp(2i7rz). 
The citation is from Iwaniec's book 

As shown in the previous section the homodyne detector behaves as a diophan- 
tine approximator of the frequency ratio v of input oscillators. The approach is 
very satisfactory in explaining frequency-locking effects in the open loop, but the 
recourse to non linear differential equations is necessary for the case of the closed 
loop. Now we attempt to develop a pure arithmetical frame to account for the 
phase-locking effects as well. This is done by normalizing the beat frequency with 




(14) 




(15) 



respect to the low pass cut-off frequency as y — instead of the reference fre- 
quency. This is suggested by the geometry of continued fraction expansions which 
resorts to new mathematical concepts such as Ford circles, the hyperbolic half-plane, 
Mobius transformations, hyperbolic Laplace equation... 

The hyperbolic half-plane is defined as the set of complex numbers z such that 
3z > 0. It is a rich mathematical object first studied by the mathematician Henri 
Poincare at the end of nineteenth century. The geometry of Ti is related to continued 
fraction expansions and it is also a natural frame to study the prime numbers. It 
will be used also for studying phase-locking effects. We put 1 . 



w UJ LF 

v + nj, v = — , < y = < 1. 



(16) 



It is clear that y > is a condition which is imposed by counting measurements. The 
second condition y < 1 results from the low pass filtering and will reveal important 
in our introduction of Ford circles below. Let us start with the continued fraction 
expansion Q which is rewritten as {eto; a±, • ■ ■ ai} = ao 

By induction po — ao, go = 1, {ao, a i} = a o 



{a 1 ;a 2 ---a i } 
apai + l _ £l 
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(17) 



Taking determinants and using a propery of continued fractions one gets ptqt-i 

I'. :<l, = i 1':' ' 
Mobius transformations 



This suggests to associate to the convergents |i the group of 



70) 



PjZ + Pi-i 

qiZ + qt-x 



PiQi-i -pi-m 



i. 



(18) 



A Mobius transform map a point z G 7i, with Ti. the upper half-plane 3z > into a 
point z' e H. Mobius transforms form a group T called the modular group. It is a 
discontinuous group that is used to tesselate Ti by copies of a fundamental domain 
T under the group action. The fundamental domain of T (or modular surface) is 
defined as T — { z £ Tt : \z\ > 1, \v\ < h}, and the family of domains {7 (J 7 ), 7 € T} 
induces a tesselation of Ti [TU] , 

We are also interested by the Ford circles they are defined by the images 
of the filtering line z = v + i under all Mobius transformations (|18|l . Rewriting 
(I18|l as qiZ 1 — Pi = —l/(qiZ + qt-i) and inserting z = v + i in that expression one 
immediately gets 



(19) 



which are circles of radius ttt centered at — 

2<3f qi 



29f 



To each 21 a Ford circle in the 



upper-half plane can be attached, which is tangent to the real axis at v = 21 . Ford 



lr The imaginary number i such that i 2 
Pi, qi and related integers. 



-1 should not be confused with the index i in integers 
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circles never intersect: they are tangent to each other if and only if they belong to 

fractions which are adjacent in the Farey sequence j < ■ ■ ■ ^ < < ^ • • • < j- 

Ford circles can be considered a complex plane view of continued fractions. 

Ford circles are also related to the wavefronts of an hyperbolic noise model: the 

Laplace equation for Ti.. To see this one first observe that Ti carries a non-Euclidean 

metric dz = (dv 2 + dy 2 ) 1 ^ 2 . This is easily shown by using ££7(2) = ^ q . z _^ q . t p and 

di{z) _ 1 cinno 1 rf- r(*0 1 _ s-yjz) 

v 

The invariance of the metric is inherited by the non-Euclidean Laplacian 



[^FFi~TP ; smct; I I = — r~ under all group actions. 



r) 2 B 2 

A = "V + v>' (20) 

We will be concerned by eigenfunctions f s (z) of (|2U[1 as our dynamical model of 
phase-locking. The complex parameter s has been introduced as a label for the 
eigenvalues. The relevant wave equation is 

(A + X)^ s (z) = 0. (21) 

The solutions of (|21() are called automorphic forms. 

They are several ways of finding eigenfunctions ^ s (z) for a given eigenvalue A. 
We will first consider elementary power law solutions, then we will show that there 
is a method to generate a lot of eigenvalues out of a fixed ^ s (z) by shifting to 
^(7(2:)), and still more by averaging over selected 7(2) in T. 

The simplest solutions of (|21fl are horizontal waves of eigenvalue A = s(s — 1) in 
the form of power laws 

*.(*) = { • (22) 

Any other wave which satisfies (|21(l can be searched by shifting z to 7(2:) so that y 
goes to 3(7(2)) with the result 

*- w - (3 *»" - krifcrF- (23) 

The wavefronts for that solution are obtained by assigning to the quantity y/(qiZ + 
<7i_i) some constant value c. As a result one gets the equation [x + qi-1/qi) 2 + (y — 
l/(2q 2 c 2 )) 2 = l/(4c 4: qf) of circles tangent to the real axis at x = —qi-i/qi and of 
radius l/(2c 2 q,^). At c = 1 such a circle is obtained from the Ford circle of index 
i by an horizontal slip from v = p%jq% to qi-i/qi- It contracts (c < 1) or expands 
(c > 1) remaining tangent to the real axis. 

The more general solution of (|21|l is obtained by summing solutions (|23|l over 7 
transformations in Y taken over the right coset r m \ T where is the subgroup 
of translations z — ► z + n. One gets 

where the summation index i means that the summation should be applied to all 
Farey fractions 



The whole solution can be decomposed into three contributions, the horizon- 
tal wave y s is incident onto the filtering line, the reflected wave S(s)y 1 ~ s is scat- 
tered with a scattering coefficient S(s) of modulus 1, whereas the remaining part 
T s (y,v) is a complex superposition of waves depending on y and the harmonics 
of exp(2i7w), but goes to zero away from the real line. More precisely one ob- 
tains S(s) = A(s)Z(s), Z(s) = , A(s) = r(1/2) r r ( ir 1/2) ' r ( s ) is thc Gamma 
function and C(s) = Y^Li ^ is the Riemann zeta function. There is much to 




Fig 2. The phase angle ft(fc) for the scattering of noise waves on the modular surface. Plain lines: 
Exact phase factor. Dotted lines: Approximation based on the quotient A(k) of two Riemann zeta 
functions at 2ik and 2ik + 1. 

say about the singularities of the Riemann zeta function. It is analytic only on 
that part of complex plane where 5Rs > 1. But in the scattering coefficient it 
appears in the extended Riemann zeta function £(s) = n~ s / 2 r(s/2)£(s) which is 
analytic over the whole complex plane. The scattering coefficient equals the quo- 
tient S(s) = £(2s — l)/£(2s). One knows from the theory of zeta function that 
s = 2 is an axis of symmetry for £(s), and according to Riemann hypothesis all non 
trivial zeros lie on that axis, with a very random distribution. An interesting case 
is thus to compute the scattering coefficient S(s) along the critical line s — \ + ik 
in which case the superposition T s (y, v) also vanishes. As a result 

S(k) oc exp(2iK(fc)) with n'ik) = dlnZ ^ &t s =-+ik. (25) 

as 2 

In that expression we removed the smooth part A(k) in the scattering coefficient 
so that its modulus is no longer 1 but its phase variability is still well described 
by the exponential factor above (see Fig. 2). The aim of this replacement is to 
get an explicit link of the counting function n'(k) to the Mangoldt function A(n) 
already accounted in (|14fl as an effective coupling coefficient in the Arnold map 
It is known that the logarithmic derivative of ((s) is — ^j~y = 5Z„>i ^T^- ^ 
is also easy to check that with B(s) — ,f , one gets = ^ n>1 ^r- with 
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b ( n j = A W(") t he modified Mangoldt function. Here Z(s) = B(2s - 1). 

It is now clear that the scattering of waves in the hyperbolic plane has much 
to share with the phase- locking model formulated in Sect. [3 One argument in 
favor of that new frame is that the average modified Mangoldt function satisfies 
B(t) — \ ~}2, n>1 b(n) = 1 + es(i), where es(i) is an error term of power spectral 
density equal to 1// 2G with G ~ 0.618 the Golden ratio [T2|. 

4. Quantum Phase-Locking 

Apparently Dirac was the first to attempt a definition of a phase operator by means 
of an operator amplitude and phase decomposition. As we have discussed, with a 
complex c-number a = Re 1 ® one obtains the phase via = aj R. Similarly, he 
sought to decompose the annihilation operator a into amplitude and phase compo- 
nents... After a brief calculation we obtain a relation indicating that the number 
operator N and phase operator $ are canonically conjugate 

[N,$\ = l. (26) 

The equation immediately leads to a number-phase uncertainty relation which is 
often seen 

6NS^>> 1/2. (27) 

However, all of the previous development founders upon closer examination. 

This is taken from |13| . a comprehensive review of the quantum phase problem. 
See also [H]. 

To approach the phase-locking problem within quantum mechanics one can start 
from the theory of the harmonic oscillator. The natural objects are the Fock states 
(the photon occupation states) |n) who live in an infinite dimensional Hilbert space. 
They are orthogonal to each other: (n\m) = 5 mn , where 5 mn is the Dirac symbol. 
The states form a complete set: J2n°=o \ n )( n \ = 1- 

The annihilation operator removes one photon from the electromagnetic field 

a\n) = y/n\n - l),n = 1,2, • • • (28) 

Similarly the creation operator cv adds one photon: a'\n) = \Jn + l\n + l),n = 
0, 1, • ■ • There is the commutation relation [a, af] = 1. The operator N — aa' has 
the meaning of the particle number operator and satisfies the eigenvalue equation 
N\n) = n\n). 

Eigenvalues of the annihilation operator are the so-called coherent states \a) 

a\a)=a\a), width |«) = e "H 2 /2 £ (29) 

71 = ^ '' 

where the eigenvalue a is a complex parameter which quantifies the intensity of 
the field. A single mode laser operated well above threshold generates a coherent 
state excitation. The electric field variation of a coherent state approaches that of a 
classical wave of stable amplitude and fixed phase. In the coordinate representation 
it is a minimal uncertainty state. 



States of well defined phase escaping the inconsistencies of Dirac's formulation 
were build by Susskind and Glogower jlfi| . They correspond to the eigenvalues of 
the exponential operator 

oo 

E = e i$ = (N + 1)-V a a = ^ \n)(n + 1|. (30) 

n=0 

Using the Hermitian conjugate operator E^ = e~ 1 *, one gets EE' = 1, E'E = 
1 — |0)(0|, i.e. the unitarity of E is spoiled by the vacuum-state projector |0)(0|. 
The Susskind-Glogower phase states satisfy the eigenvalue equation E\^>) = e l ^\^>); 
they are given as 

oo 

|*> = Xy^|n>. (31) 

?i=0 

Like the coherent states the phase states are non orthogonal and they form an 
overcomplete basis which solves the identity operator: ^ $*dip\E){E\ — 1. The 
operator cos$ = ^(E + E^) is used in the theory of Cooper pair box with a very 
thin junction when the junction energy Ej cos $ is higher than the electrostatic 
energy [T%| . 

Further progress in the definition of phase operator was obtained by Pegg and 
Barnett. The phase states are now defined from the discrete Fourier transform (or 
more precisely the quantum Fourier transform since the superposition is on Fock 
states not on real numbers) 

1 9-1 

\9 P ) = — ^exp(2wr-ra)|n). (32) 
n=o q 

The states are eigenstates of the Hermitian phase operator 

9-1 

9« = I>pIWpI, (33) 

with 6 p = 9q + 2irp/q and #o is a reference angle. It is implicit in the definition 
1(55)1 that the Hilbert space is of finite dimension q. The states \9 P ) form an or- 
thonormal set and in addition the projector over the subspace of phase states is 
Sp=o \&p){8p\ = Igj where l q is the unitary operator. Given a state \F) one can 
write a probability distribution |(0p|-F)| 2 which may be used to compute various 
moments, e.g. expectation values, variances. The key element of the formalism is 
that first the calculations are done in the subspace of dimension q, then the limit 
q — * oo is taken |17j. 

We are now in position to define a quantum phase-locking operator. Our view- 
point has much to share with the classical phase-locking problem as soon as one 
reinterpret the fraction ^ in i[32J) as arising from the resonant interaction between 
two oscillators and the dimension q as a number which defines the resolution of the 
experiment. From now we emphasize such phase states \9' p ) which satisfy phase- 
locking properties and we impose the coprimality condition 



(34) 
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The quantum phase-locking operator is defined as 

e ^ = j29 p \e' p ){e' p \, (35) 

p 

with Bp = 27r^ and the notation p means summation from to q — 1 with (p, q) = 1. 
Using (j22J) and (PU in one obtains 

Q^ = lj2c q (n-l)\n){l\, (36) 

9 n,l 

where the range of values of n, I is from to <j){q), and 4>{q) is the Euler totient 
function. The coefficients in front of the outer products \n)(l\ are the so-called 
Ramanujan sums 

c q( n ) = Vexp(2i7r-7i) = ^i)<Hg) with = q/( q n )_ (37) 

In the above equation //(g) is the Mobius function, which is if the prime number 
decomposition of q contains a square, 1 if q = 1 and (— l) fe if q is the product of 
k distinct primes. Ramanujan sums are relative integers which are quasi-periodic 
versus n with quasi-period <p(q), and aperiodic versus q with a type of variabil- 
ity imposed by the Mobius function. Ramanujan sums have been used for signal 
processing of low frequency noise ^]. With the Ramanujan sum expansion the 
modified Mangoldt function introduced at the end of Sect. El is the dual of Mobius 
function 

H^%) = E^W (38) 

This illustrates that many "interesting" arithmetical functions carry the structure of 
prime numbers. We already mentioned the relation d In ((s)/ds = J2 n >i ^J^i 
there is also the relation 1/C(s) = Y^ n >i ^« • There is a well known formulation of 
Riemann hypothesis from the summatory Mobius function 52„_i A*( n ) = 0(t 1 / 2+e ), 
whatever e pQ. 

Given a state (3 one can calculate the expectation value of the quantum phase- 
locking operator as 

p 

If one uses the finite form of Susskind-Glogower phase states Ij31|) and a real pa- 
rameter /3 

1 9-1 

|/?) = — 5>xp(m/?)|n), (40) 

V ? n=0 

the expectation value of the locked phase becomes 



5 10 15 20 25 30 35 40 45 50 
dimension q 



Fig 3. Oscillations in the expectation value 1411 of the locked phase at /3 = 1 (dotted lines) and 
their squeezing at (3 = (plain lines). The brokenhearted line which touches the horizontal axis 
is 7rA(g) / In q 

For (3 = 1 it is found that (Q l ° ck ) has the more pronounced peaks are at such values 
of q which are powers of a prime number. It can be approximated by the normalized 
Mangoldt function irA(q)/ \nq as shown on Fig. 3. For f3 = the expectation value 
of (8q° ck ) is much lower. The parameter (3 can be used to minimize the phase 
uncertainty well below the classical value [3] . 

A remarkable model of the quantum phase-locking effect and its relation to 
prime number theory has been constructed by Bost and Connes |2()| . Instead of an 
ad-hoc quantum phase operator as ( 17 1) or J2EJ), it is based on the formulation of 
a dynamical system and its associated quantum statistics. The dynamical system is 
first defined by an Hamiltonian operator Hq with eigenvalues equal to the logarithms 
of integers 

H \n) = hm|n). (42) 

Using the relations exp(— (3Ho)\n) — exp(— (3lnn)\n) = n _/3 |n), it follows that the 
partition function of the model at the inverse temperature (3 is 

oo 

Trace(exp(-/3H )) = £ n^ = (0), (43) 

n=l 

where Cift) is the Riemann zeta function already met in Sec. [21 in the definition of 
the hyperbolic scattering coefficient (|25(l . In quantum statistical mechanics, given 
an observable Hermitian operator M one has the Hamiltonian evolution at{M) 
versus time t 

a t (M) = e itHo Me~ itHo , (44) 
and the expectation value of M is the Gibbs state 

Gibbs(M) = Trace(Mexp(-/3H ))/Trace(exp(-/3H ). In Bost and Connes ap- 
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Fig 4. Phase expectation value 1471 in Bost and Connes model at the inverse temperature /3 = 3 
(plain lines) in comparison to the function fi(q) / <f>(q) (dotted lines) 
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Fig 5. Phase expectation value 1471 in Bost and Connes model at the inverse temperature j3 = 
e = 0.1 (plain lines) in comparison to the function —A(q)e/q (dotted lines). 



preach the observables belong to an algebra of operators fi q and e q which are 
defined by their action on the occupation numbers \n) as 

fi q \n) = \qn), (45) 
4 P V> = exp(2Hp)|n). (46) 

The first operator /i q acts as a shift in the space of number states; the second one is 
such that its action encodes the individuals in the quantum Fourier transform (|32[) . 
Like in the quantum phase-locking operator one uses the coprimality condition l|34() 
to distinguish in (|46|l the primitive roots of unity exp(2iirp/ q) , (p,q) = 1. One can 
show that there is a hidden symmetry group which is used to label the elements 
of the algebra 2 . Using the action of the group, the Gibbs state is replaced by the 
so-called Kubo-Martin-Schwinger (or KMS) state. The system exhibits a phase 
transition with spontaneous symmetry breaking at the inverse temperature (3=1 
which corresponds to the unique pole of the Riemann zeta function £(/3). At low 
temperature /3 > 1 one gets, after tricky calculations, the expectation value of the 
phase operator which replaces (|41J) in the following form [20] 

KMS(e(P>)=q^ JJ (47) 

p prime 

The KMS state is represented for two limiting cases, the low temperature limit 
/9 > 1 (Fig. 4) and the critical case j3 = 1 + e (Fig. 5), with e ~ 0. In these 
limits one has respectively KMS^»i(q) = and KMSi +e ~ — ^si£. T n tfie 

low temperature limit the spectrum l|38|l corresponding to the Ramanujan sum 
expansion of the modified Mangoldt function b(n) is recovered. Close to the critical 
point /? = 1 + e the oscillations are proportional to A(q) ~ b(q) and are of very 
small amplitude due to the squeezing coefficient e. A comparable squeezing effect 
was already observed in the expectation value (|39|1 of the quantum phase operator 
(see Fig. 3). After the phenomenological model (|15|l and the hyperbolic model Q25JI. 
the Bost and Connes cyclotomic model also points to the Mangoldt function as a 
source of low frequency fluctuations. In the last case the model is associated to 
the spontaneous symmetry breaking and the squeezing of phase oscillations at the 
critical KMS state. 

As a conclusion phase-locking evokes Einstein "spooky" action at a distance. It 
is a very important concept in phase sensitive communication circuits. Quantum 
phase- locking may become synonymous of phase entanglement, be used to such tasks 
as remote synchronization, quantum noise reduction and as a quantum resource for 
the circuits of quantum communications. 
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2 This is the Galois group W = Gal(Q cyci / Q) of the cyclotomic extension on the field of rational 
numbers Q. 
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